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We present a method for obtaining the quasi-exact solutions of the Rabi 
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1 Introduction: formal aspects 

A new class of potentials which are intermediate to exactly solvable (ES) ones and 
non solvable ones are called quasi-exactly solvable (QES) problems [H [2j [3l H]. It 
is well known that a part of their spectrum can be determined algebraically but not 
whole spectrum. One of an alternative mathematical language for describing the QES 
systems is a remarkable set of orthogonal polynomials introduced by Bender and Dunne 
[5]. These polynomials satisfy the three-term recursion relation and, as a consequence, 
form an orthogonal set with respect to some weight function p(S) [HI U\- The exact 
solution of the Schrodinger equation takes the form: 
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in which V n (£) denotes certain polynomials of £ and satisfy [H [9] : 

V n (£ k ) = 0; n>A+l, k = 0, 1, 2, ■ ■ ■ , A. (2) 
where the parameter A is a positive integer value. It follows that V n (£) are orthogonal: 

V n {£) p(£)V m (£)d£ = PnP m 5 nm , (3) 



where p n represent the norms of the orthogonal polynomials V n (£)■ It is possible 
to determine the norms p n (or squared norms 7 n = p 2 n ) by multiplying the recursion 
relation by £ n ~ x p (£) and to integrate with respect to £. 

However, the search of the orthogonal polynomials V n (£) can be generated by 
means of the asymptotic iteration method (AIM) [10J. This procedure was first intro- 
duced as an approximation in order to deduce both energy eigenvalues and eigenfunc- 
tions, using some computer-algebra systems. Due to the simplicity of the procedure, 
various aspects of the method have been employed, successfully, to obtain both ES 
and QES problems [JUl HH [12], [131 EE]- The systematic procedure of the AIM begins 
by rewriting a second-order linear differential equation in the following form: 

X" (x) = r (x) x (x) + s (x) x (x) , (4) 

with ro (x) and sq (x) are many times differentiate functions. In order to find the 
solution of Eq.(4), we rely on the symmetric shape of the right-hand side. Thus for 
(n + l) th and (n + 2) th derivative of Eq.(4), n — 1, 2, • • ■ , we get: 

r n (x) = r' n _ x (x) + s n _i (x) + r (x) r n _i (x) , (5.1) 
s n (x) = 4-i ( x ) + s o ( x ) r n-\ (x) , (5.2) 

where for sufficiently large iteration number n, the " asymptotic" behavior of the pro- 
cedure can be applied as: 

r n _x (x, £) _ r n (x, £) _ ^ ^ 
s n _i{x 1 £) s n (x,£) 

and which allows us, on the one hand, to calculate the energy eigenvalues by iterations 
and, on the other hand, to express the so-called polynomial solutions of Eq.(4). 

Many interesting models are obtained by combining two fundamental models of 
quantum mechanics; namely the two-level system and the harmonic oscillator [T5] . 
The Rabi Hamiltonian is one of the more studied ones [TBI [T7] , and plays an important 
role in many areas of physics from condensed matter and biophysics to quantum optics. 
Due to the elegance of the model various aspects of the Hamiltonian have been studied 
analytically and numerically [HI [191 EQl EU [221 E3]. The complete Hamiltonian of a 
such system is of the form: 

HQ 

H = hula) a + -—-co + A (cr + + cr_) (a' + a) , (7) 
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where the parameter ou is the boson field frequency, HVt is the atomic level separation 
and A is the atom-field coupling constant. Here a) (resp. a) is a creation (resp. an- 
nihilation) operator for the field mode and the a are the usual Pauli spin matrices 
o x = a + + cr_, o y = —i (cr+ — a ) and a z = a . Note that the zero energy level was 
taken halfway between the two atomic levels, so that the unperturbed atomic ener- 
gies are ±4p. The four terms appearing in the interaction part of the Hamiltonian 
(7) have the following interpretation: a + a (resp. a^a^) represents an absorbed (resp. 
emitted) photon and an excited (resp. de-excited) atom from state |1) to |2), while 
a+al (resp. cr_a) stands for one photon which is emitted (resp. absorbed) and an 
excited (resp. de-excited) atom. 

In this present paper, we propose an alternative approach to study the Rabi 
Hamiltonian (7) based to the AIM procedure, to obtain both the energy eigenvalues 
S n , the eigenfunctions x (£j £) an d a set of orthogonal polynomials V n (£) explicitly 
constructed from the above procedure and known to be orthogonal polynomials of 
Bender-Dunne of the genus one. In particular, we show firstly that these orthogonal 
polynomials generated from the three-term recursion relation and the coefficients of 
the polynomial solutions deduced starting from AIM procedure are equivalent, and 
secondly that the orthogonal polynomials V n (£) are prone with a nonpositive definite 
norm, unless the constraint n > E + A 2 + 1 is checked. 

The paper proceeds as the following. In the next section we introduce a unitary 
transformation with an aim of making the Rabi Hamiltonian simpler. By applying 
AIM procedure, we deduce the coefficients of the polynomial solutions with the first 
five iterations. In section 3 we clarify, through an explicit construction, the relation 
between the calculated eigenfunctions and the orthogonal polynomials leading to the 
three-term recursion relation and which enables us to deduce the orthogonal polyno- 
mials V n {£)■ We show, in section 4, that the orthogonal polynomials are associated 
with a nonpositive norm 7 n and the last section is devoted to our final conclusion. 



2 Solvability of the Rabi Hamiltonian 

In Fock-Bargmann space, the system (7) is transformed into the third-order differen- 
tial equation which is not easy to solve. Therefore, the first step to solve the system 
(7) is to reduce the order of the differential equation using some transformations. 
The structure of the Rabi Hamiltonian is more easily seen by applying the unitary 
transformation: 

H = U ] HU 

= a t a + /3(cr + + cr_) + Ao- (a t + a) , (8) 
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where U = — 4g (a + + a_ + a ) . Here, fi and u; are set to unity, and Q = 2/3. By 
performing a transition to a Fock-Bargmann space — >• x and a — >■ -p, and writing 



the stationary Schrodinger equation for the two-component wave- function (Jig), we 
obtain a system of two first order linear differential equations for the functions if) 1 (x) 
and ip 2 (x) [HI [20]: 

dipt (x) £ - Ax /S . . , 

ax x + A x + A 

<M 9 (x) /9 , . x 8 + Ax , . s . 

^jr 1 = — ix + — r ^ 2 * , 9.2 

ax x — A x — A 

where £ are energy eigenvalues of %. The general solution of Eqs.(9) can be obtained 
by transforming them into a second-order form: 

(x 2 - A 2 ) ifj'l (x) - {{£ - Ax - 1) (x - A) + (x + A) {£ + Ax)] $ x (x) 

- [/3 2 - £ 2 + A 2 x 2 - A (x - A)] i) x (x) = 0, (10) 

and substituting 

x = A(2£-l), ^( X ) = e- 2X \(0, (11) 
we obtain a second-order differential equation [i~9] : 

C(l-0x"(0+ [A 2 (4^ 2 - 2^ - 1) + ^ (2^ - 1) - e + 1] 

+ [A 4 (3 - 40 + 2£A 2 (1 - 20 + /? 2 - £ 2 ] x (£) = 0. (12) 

Now this equation is similar to Eq. (4) and then it is in suitable form for application 
of AIM procedure; the initial ro (£) and so (£) functions are given: 

ro(0 = A^-2 gz ^ (131) 
^ _ A 4 (3 - 4fl + 2£A 2 (1 - 2fl + I 2 - £ 2 

S 0UJ - " U^.^J 

We may calculate r n (0 and s n (0 using a sequence of Eqs.(5) and the calculated 
energy eigenvalues £ n by the mean of Eq.(6) should be independent of the choice of 
the variable £. 

We report below the factored form of the first five iterations of Eq.(12) which lead to 
the coefficients of the polynomial solutions. By convention n represents the iteration 
number, while d refers to the degree of the polynomial solutions. The coefficients 
y n (£; A, f$) are polynomials in the energy variable £ n with an even higher-degree 
2n + 2. 

Note that the energy eigenvalue index is the same as the iteration number. 
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n 


= 1, 


d = 


2, 


Cl2 


= 16A 4 (-l+£ + A 2 ) (£ + A 2 ) 






d = 


1, 


Cu 


= 8A 2 (-l+£ + X 2 )y (£;X,/3) 






d = 


o, 


Cio 


= yi(e;\j) 


n 


= 2, 


d = 


3, 


C23 


= 64A 6 (-2 + £ + A 2 ) (-1 + £ + A 2 ) (£ + A 2 ) 






d = 


2, 


C22 


= 48A 4 (-2 + £ + A 2 ) (-1 + £ + A 2 ) y (£; A, 0) 






d = 


1, 


C21 


= 12A 2 (-2 + £ + A 2 ) ^ (5; A, (3) 






d = 


o, 


C20 


= ^ 2 (^;A,/3) 


n 


= 3, 


d = 


4, 


C34 


= 256A 8 (-3 + £ + A 2 ) (-2 + £ + A 2 ) (-1 + £ + A 2 ) (£ + A 2 ) 






d = 


3, 


C33 


= 256A 6 (-3 + £ + A 2 ) (-2 + £ + A 2 ) (-1 + 5 + A 2 ) y (5; A, f3) 






d = 


2, 


C32 


= 96A 4 (-3 + £ + A 2 ) (-2 + £ + A 2 ) ^ (5; A, (3) 






d = 


1, 


C31 


= 16A 2 (-3 + £ + A 2 ) y 2 (£; A, /3) 






d = 


o, 


C30 


= y 3 (£;A,/3) 


n 


= 4, 


d = 


5, 


C45 


= 1024A 10 (-4 + £ + A 2 ) (-3 + £ + A 2 ) (-2 + £ + A 2 ) 
(-l + £ + A 2 ) (£ + A 2 ) 






d = 


4, 


C44 


= 1280A 8 (-4 + £ + A 2 ) (-3 + £ + A 2 ) (-2 + £ + A 2 ) 

(-i + £ + A 2 );yo(£;A,/3) 






d = 


3, 


C43 


= 640A 6 (-4 + £ + A 2 ) (-3 + £ + A 2 ) (-2 + £ + A 2 ) ^ (5; A, (3) 
= 160A 4 (-4 + £ + X) (-3 + £ + A 2 ) y 2 (£; A, /3) 






d = 


2, 


C42 






d = 


1, 


C41 


= 20A 2 {-4 + £ + X 2 )y 3 (£;\,f3) 






d = 


o, 


C40 


= y,(£;X,P) 


n 


= 5, 


d = 


6, 


^56 


= 4096A 12 (-5 + S + A 2 ) (-4 + £ + A 2 ) (-3 + £ + A 2 ) 
(-2 + £ + A 2 ) (-l+£ + A 2 ) (£ + A 2 ) 






d = 


5, 




= 1280A 10 (-5 + £ + A 2 ) (-4 + £ + A 2 ) (-3 + £ + A 2 ) 
(-2 + £ + A 2 ) (-1 + £ + A 2 ) 3^o (£; A, 0) 






d = 


4, 


C54 


= 3840A 8 (-5 + £ + A 2 ) (-4 + £ + A 2 ) (-3 + £ + A 2 ) 
(-2 + £ + A 2 )^(£;A,/3) 






d = 


3, 


C53 


= 1280A 6 (-5 + £ + A 2 ) (-4 + £ + A 2 ) (-3 + £ + A 2 ) y 2 {£■ A, /?) 






d = 


2, 


C52 


= 240A 4 (-5 + £ + A 2 ) (-4 + £ + A 2 ) y 3 (5; A, 0) 






d = 


1, 


C51 


= 24A 2 (-5 + £ + A 2 ) 3^4 (5; A, (3) 






d = 


o, 


^50 





These polynomials have some interesting properties. For example, a polynomial set 
y m (£] X, (3) can be generated successively for each a new iteration with the constraint 
n = d + m, where < m < n. However if n — d (i.e. m — 0), then the associated 
polynomial solution is: 



y (£;X,/3) = £ 2 -2X 2 £- (3 2 -ZX\ (14) 

The factor associated to the coefficient C n ,n+i is (4A 2 ) n+1 and C„o = y n (£; X, /3), 
with n — 1,2,.... Furthermore, we will see that the polynomials y n (£; A, (3) lead to 
the Bender-Dunne orthogonal polynomials as reviewed in the next section. 
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3 Bender— Dunne orthogonal polynomials 



The search for a power series solution of Eq.(12), x (0 = J2^=o Xn (£)> can be explicitly 
generated by means of the series expansion leading to the three-term recursion relation. 
Using the identities: 

oo 

x(0 = ? 9 $>^ n ' ( m ) 

n=0 

oo 

n=0 

oo 

X"(0 = ^^(n + g)(n + g-l)a n r- 2 , (15.3) 

ra=0 

where the exponent q and the coefficients a n are still undetermined, and by substituting 
them into Eq.(12), we have: 

oo oo 

AX 2 (n + q - X 2 - S) a n C +q+l + ]T [3A 4 - S 2 + (3 2 + 2SX 2 

n=0 n=0 

oo 

+ (n + ?) (25 - A 2 - n - q)] a n C +q + £ (n + q) {n + q - A 2 ) a^^" 1 = 0, (16) 

n=0 

where the coefficients of the left-hand side of Eq.(16) must vanish individually. The 
lowest power of £ appearing in Eq.(16) is £ 9_1 , for n = in the last summation. The 
requirement that the coefficients vanish yields the indicia! equation q [q — A 2 ) a = 0; 
hence, we must require either q = or q = A 2 , and a ^ 0. Note that if the first 
summation vanishes, the energy eigenvalues satisfy the constraint: 

£ n , q (X) = n + q-X 2 . (17) 

Therefore, if we replace n — > j in the second summation and (n — 1) — > j in the last 
(they are independent summations), this results in the two-term recursion relation: 

%+i - - j—, ; i w ■ ; n T2V a i- ( 18 ) 

(j + q + 1) (j + q + 1 - A j 

Substituting j by n and choosing the indicial equation root q — 0, the mathematical 
induction leads us to write: 

(_l)™/3 2 " T(l-A 2 ) 

ra! T(n + 1-A 2 ) K J 

with a = 1. The solution of Eq.(12), following Eq.(l) 

oo 

X&S) = £Xn(0^(£) 
n=0 

_ - (-iyg; r(i-A') 

2^ n! Wn + i-A 2 ) " l K • ( 2U » 

n=0 V / 
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can be considered as the generating function for the polynomials V n (£). However, for 
positive integer values A, the series expansion in Eq.(20) are truncated when £ is a zero 
of V\ (£). We know that r (n + 1 — A 2 ) has simple poles at n = A 2 — 1, A 2 — 2, A 2 — 3, . . ., 
[21] and taking into account the constraint (17) with q = 0, leading to the upper limit 
for n, i.e. n max = A = £ + A 2 > 0. Indeed, this can be observed in the coefficients C nc i 
of the polynomial solutions like it was repoted above. Therefore, the corresponding 
exact eigenfunctions are given: 

rc=0 V / 

By substituting Eq.(21) into Eq.(12) leads to the following three-term recursion 
relation for V n {£): 

4A 2 (n-e- A 2 ) V n+1 {£) + [(n-£- A 2 ) {£ - 3A 2 - n) - f3 2 ] V n {£) 

+ 4\ 2 (n-£-X 2 )V n - 1 (£) = 0, (22) 

with the normalization condition Vq {£) = 1. The polynomial V n (£) vanishes for 
n > £ + A 2 + 1 as occurred in the third term of Eq.(22). The recursion relation (22) 
generates a set of polynomials, where the first five of them are reported below in table 
1. 



n 


£(X) 




1 


1 - 


A 2 


4A 2 + (/? 2 - 1) 


2 


2 - 


A 2 


32A 4 + 4 (3/3 2 - 8) A 2 + (/f - l) (/3 2 - 4) 


3 


3- 


A 2 


384A 6 + 16 (11/3 2 - 54) A 4 + 8 (3/3 4 - 29/3 2 + 54) A 2 
+ (/3 2 - 1) (/3 2 - 4) (/? 2 - 9) 


4 


4- 


A 2 


6144A 8 + 128 (25/3 2 - 192) A 6 + 16 (35/3 4 - 542/3 2 + 1728) A 4 + 8 (5/3 6 
-115/3 4 + 722/3 2 - 1152) A 2 + (/3 2 - l) (/3 2 - 4) (/3 2 - 9) (/3 2 - 16) 


5 


5- 


A 2 


122880A 10 + 512 (137/? 2 - 1500) A 8 + 64 (225/3 4 - 50 36/3 2 + 24 00) A 6 
+16 (85^ 6 - 28 6 7/3 4 + 275 1 8/? 2 - 72000) A 4 + 4 (l5/? 8 - 670/5 6 + 9551/3 4 
-49216/3 2 + 72000) A 2 + (/3 2 - l) (,5 2 - 4) (f3 2 - 9) (/3 2 - 16) (/3 2 - 25) 



Table 1: Bender-Dunne orthogonal polynomials for Rabi Hamiltonian and the associ- 
ated energy eigenvalues 



These orthogonal polynomials are exactly similar to those obtained by the Juddian 
isolated method and agree with the first roots of the Kus series pT8l fT9] . On the other 
word, the system has energy £ n = n — A 2 , n = 1,2,..., only if the atomic-level 
separation 2(3 and the boson-atom field coupling A obey condition V n (A) = 0. Note 
also that these orthogonal polynomials can be deduced starting from the coefficients 
C n d of the polynomial solutions y n (£; A, 0) with the same associated energy eigenvalues 
(17). 
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4 Norms of orthogonal polynomials 



As the polynomials V n (£) are orthogonal of a discrete variable £ n , it is possible to 
determine their squared norms. The procedure consists to apply Eq.(3); i.e. multi- 
plying the recursion relation (22) by £ n ~ l p(£), where p(£) is the weight function, 
and to integrate with respect to £ using the fact that both V n (£) and £ k , k < n, are 
orthogonal. We obtain: 

, nfl — A 2 , . 

7n = 2A 2 — — ^n-i, (23) 
n + A 

with 7 n = p 2 and 7 = 1. It is obvious that by mathematical induction, we get: 

! , k-X 2 
k + X 2 



7n = Il- V -~ 

k=l 



9 (1-A 2 ) 

= 2"A 2 " | i + A2 j" , (24) 

where (a) = a (a + 1) (a + 2) • • • (a + n — 1) are the often-used Pochhammer symbol 
[21], with (a) = 1. Using the well-known identity r (z) T (1 — z) = [21], we 

have: 

9 n X 2n sinTrA^ , 2 , , 2 v r (w + 1 - A 2 ) 
7„ = 2A -^r(A)r(l + A) r( ^ + i + A2) . (25) 

This identity reveals an important result that the orthogonal polynomials are 
associated with a nonpositive definite norm represented by r (n + 1 + A 2 ) . How- 
ever, it shows in the T-function of the numerator that 7 n has in the complex plan 
poles which are given at the negative integers iV = 0,-1, —2, . . ., (or equivalently: 
n = A 2 — 1, A 2 — 2, . . .), leading to the upper limit A = £ + A 2 > as conjectured 
above. Consequently, 7 n vanishes for each n > A + 1. 



5 Conclusion 

In this paper, we have shown that there exists another approach of introducing the 
Bender-Dunne polynomials which characterizes the quasi-exact solvable problems. 
The model given here is related to the AIM procedure and we have taken a new look 
at the solution of the Rabi Hamiltonian in order to generate its energy eigenvalues as 
well as the eigenfunctions. 

The main aim of this article was to consider the possibility of generating the 
Bender-Dunne polynomials V n (£) via two different but equivalent ways. The first 
way, considered as the standard method, consists in applying the eigenfunctions to 
the three-term recursion relation, while the second enables us to generate them di- 
rectly starting from the coefficients C n a of the polynomial solutions y n (£; A, f3) deduced 
by AIM procedure. Consequently, it proves that V n (£) are similar to 34 (£] A, j3) if 
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the system has the associated energy eigenvalues £ n = n — A 2 ,n=l,2,... We have 
also shown that the orthogonal polynomials arising from a quasi-exactly solvable Rabi 
model is associated with a nonpositive definite norm whose the iteration number n 
exceeds a critical value A mentioned above, i.e. n > A + 1. Consequently, the quasi- 
exact energy eigenvalues of the Rabi Hamiltonian are the zeros of the polynomials 
Vk{£). 

For other applications, we believe that similar results can be obtained via the 
approach worked out here to the various atomic systems such as Jaynes-Cummings 
and 8 <E> e Jahn-Teller Hamiltonians. 
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